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It is shown that transverse momentum dependent azimuthal spin asymmetries suffer from 
suppression due to Sudakov factors, in the region where the transverse momentum is much 
smaller than the large energy scale Q 2 . The size and Q 2 dependence of this suppression are 
studied numerically for two such asymmetries, both arising due to the Collins effect. General 
features are discussed of how the fall-off with Q 2 is affected by the nonperturbative Sudakov 
factor and by the transverse momentum weights and angular dependences that appear in 
different asymmetries. For a subset of asymmetries the asymptotic Q 2 behavior is calculated 
analytically, providing an upper bound for the decrease with energy of other asymmetries. 
The effect of Sudakov factors on the transverse momentum distributions is found to be very 
significant already at present-day collider energies. Therefore, it is essential to take into account 
Sudakov factors in transverse momentum dependent azimuthal spin asymmetries. 

13.65.+i; 13.88.+e 



I. INTRODUCTION 



In this article we will study the effects of Sudakov factors in transverse momentum dependent azimuthal spin 
asymmetries, like the Collins effect asymmetry of Ref. [jij. We will demonstrate explicitly that such asymmetries 
suffer from suppression due to these Sudakov factors, in the region where the transverse momentum is much 
smaller than the large energy scale Q 2 . This Sudakov suppression stems from soft gluon radiation and increases 
with energy. It implies that tree level estimates of such asymmetries tend to overestimate magnitudes and 
increasingly so with rising energy. In this paper, the Q 2 dependence of the suppression due to Sudakov factors 
will be studied numerically in two examples which are relevant for present-day experimental studies. In addition, 
the asymptotic Q 2 dependence of an important subset of asymmetries is calculated analytically, providing an 
upper bound for the decrease with energy of other asymmetries. 

The first example we consider is a Collins effect driven cos(20) asymmetry in electron-positron annihilation 
into two almost back-to-back pions ^|3|, which in principle can be determined from existing LEP data. This 
asymmetry allows for a determination of the Collins effect fragmentation function, which in turn would be 
useful for the extraction of the transversity distribution function h\. The latter can be done via the single spin 
asymmetry in semi-inclusive Deep Inelastic Scattering (DIS) which was originally proposed by Collins |lj] -our 
second example-. In Ref. p] Collins already remarks that the Sudakov factors will have the effect of diluting 
this single spin asymmetry due to broadening of the transverse momentum distribution by soft gluon emission. 
We will study this effect in a quantitative way to gain insight into the parametric dependences of such Sudakov 
suppression. 

For this we will follow the recent analysis 0] of a hclicity non-flip double transverse spin asymmetry in 
vector boson production in hadron-hadron scattering. That asymmetry is for instance relevant for the polarized 
proton-proton collisions to be performed at BNL-RHIC. In Ref. Jij the effect of the Sudakov factors compared 
to the tree level asymmetry expression was numerically estimated. It was shown that the inclusion of Sudakov 
factors cause suppression by at least an order of magnitude compared to the tree level result at scales around 
M\y or Mz- Moreover, the suppression increases with energy approximately as Q 6 (in the studied range of 
roughly 10— 100 GeV). The conclusion was that the Sudakov suppression together with a kinematic suppression 
(due to explicit lightcone momentum fractions that appear in the prefactor) imply that the asymmetry will be 
negligible for Z or W production and is interesting only at much lower energies. 

Similarly, Sudakov suppression will turn out to be an important issue for the actual determination of the 
above mentioned cos(2</>) asymmetry from LEP data (or in general at high energy e + e~ colliders). Due to the 



lack of knowledge of the nonperturbative Sudakov factor in the case of electron-positron annihilation into almost 
back-to-back hadrons, solid numerical conclusions about the size and Q 2 dependence of the suppression cannot 
be drawn. Nevertheless, one can draw conclusions about what determines the Q 2 dependence of the transverse 
momentum distribution of the asymmetries and about the size of the suppression for generic nonperturbative 
Sudakov factors. For that purpose comparisons to tree level are most instructive, since the evolution of the often 
unknown distribution and fragmentation functions themselves does not affect the relative Sudakov suppression 
compared to tree level. 

The two examples to be investigated here -the above mentioned cos(20) asymmetry and the original Collins 
sin(0) asymmetry in semi-inclusive DIS-, have different transverse momentum weights and angular dependences, 
producing different Sudakov suppression effects. General conclusions about how these properties affect the Q 2 
dependence can be drawn. In general, a larger power of transverse momentum in the weight of an asymmetry 
implies larger suppression. Both the numerical calculations at realistic collider energies and the asymptotic 
behavior of asymmetries exhibit this feature. 

Before going into the details of the specific examples, we will first give an overview of the issues regarding 
factorization theorems for transverse momentum dependent cross sections, in the region where the transverse 
momentum is very small compared to the hard scale in the process ^-Q). This will provide the theoretical 
justification of the factorized asymmetry expressions that will be derived. 

The outline of this paper is as follows. In section II we will discuss the essentials of factorization theorems for 
the Drell-Yan process (relevant for the double transverse spin asymmetry studied in Ref. J0), which are in fact 
completely analogous to the case of back-to-back jets in electron-positron annihilation [pE] and semi-inclusive 
DIS (SIDIS). Following the approach of Ref. [Q, we will then study the cos(20) asymmetry in electron-positron 
annihilation (section III) and the Collins sin(</>) asymmetry in semi- inclusive DIS (section IV). In section V we 
will discuss the asymptotic behavior of the Sudakov suppression for certain subsets of asymmetries. 



II. FACTORIZATION AND TRANSVERSE MOMENTUM 

For defmiteness, we will consider here the Drell-Yan process (Hi + H 2 — > 1 + 1+ X) and discuss what is known 
about the factorization of this process. First of all, it is well-known that a local operator product expansion 
(OPE) cannot be applied, since the process receives relevant contributions off the light cone i.e. from matrix 
elements like 

^(Pi,P 2 |J M (OI^)(^l^(0)|Pi,P 2 ) with £ 2 ,4 0. (1) 
x 

Therefore, one usually considers nonlocal operators and what is called a "working redefinition of twist" M, 
which means the lowest value of t in l/Q* -2 at which a function can contribute to the cross section (Q is the 
invariant mass of the lepton pair). The nonlocal operators appear in the leading twist factorization theorem 

, n , dg , =Y.f dx I dx$ a (x)H ab (x,x-,Q 2 )W(x), (2) 
dildxidx 2 ~^Jx! Jx 2 

where a, b are flavor indices and H ab (x, x;Q 2 ) is the hard (partonic) part of the scattering, which is a function 
of the hard scale Q 2 and the lightcone momentum fractions x,x only. The correlation functions <E> (<E> a = <E> a ) 
describe the soft (nonperturbative) physics and are given by nonlocal operator matrix elements 

! ^ Xx (P,S\^(0)C + [0,Xn^]M^-)\P,S), (3) 



where the path-ordered exponential, 



£+[0, An_] = 7>cxp (-ig J di]A + (r)n-)j , 



(4) 



renders the matrix element color gauge invariant. The hadronic state \P, S) is determined by the hadron 
momentum P and spin vector S. Also, we have suppressed the factorization scale dependence in Eq. (|^). 

Secondly, if one considers the cross section differential in the transverse momentum q T of the lepton pair, 
then for the case of q\ = Q\ <gC Q 2 collinear factorization, like in Eq. (|^), does not hold. Instead, the leading 
twist factorization takes the form P-ffpjJlOl] 
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da 



d£ldx\dx2d 2 qT 



E 



dx 



dx / d krd pr 



d 2 b 

(2tt)* 



-ib-(p T +k T -q T ) 



x$ Q (x, k T ) H ab (x, x; Q) e ~ s ^ <S>\x,p T ) + Y(x 1 ,x 2 , Q, Q T ). 



(5) 



The correlation function <j> a (a;, k T ) is now also a function of transverse momentum, the factor eT s ( b > < 2) is the 
so-called Sudakov form factor and the factor Y{x\, x%, Q, Qt) becomes important only when Qt ~ Q- Here we 
use a different notation than in e.g. Refs. [^-[7||9|,[l0t] , but to make a connection with for instance the notation 
of Ref. m we note that $>(x, k T ) corresponds to the Fourier transform of VA/i( x A> V) P aI^ an d the e~ s terms 
correspond to each other. 

The factor Y{x\, X2, Q, Qt) is present to retrieve collinear factorization when Qt is a hard scale itself, i.e. 
Qt ~ Q- We will neglect the Y term from here on, since we will only be interested in the region of Q T <C Q 2 . We 
will now comment on the other terms in the above factorization theorem and explain the additional restrictions 
we will impose in the limit Q T <§C Q 2 . 

The factor e~ s ^ b '^ is the so-called Sudakov form factor, which arises due to exponentiation of soft gluon 
contributions. This is in contrast to inclusive cross sections, like Eq. (^]), in which there is a cancellation of soft 
gluon contributions. At values b 2 = b 2 <C 1/A 2 , the Sudakov form factor is perturbatively calculable and is of 
the form 



5(6, Q) 



d^ 2 



bl/b 2 V- 



A(a„{fi)) In 
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M 2 



B{a s {n)) 



(6) 



where bo = 2exp(— 7#) w 1.123 (we choose the usual constants C\ = bo, C2 — 1). One can expand the functions 
A and B in a s and the first few coefficients are known, e.g. [ fj"l]Jl2]| . In order to obtain a first estimate of the 
effect of including the Sudakov factor we will consider only the leading contribution, i.e. take into account only 
the first term in the expansion of A: AW = C f /tt, which is the same for unpolarized as well as polarized 
scattering. This leads to the expression jl3| 



S(b,Q) 



16 



33- 2n f 



log 



b 2 Q 2 



]«>;■ ( ^_ ) | ( ,o 



log (b 2 Q 2 /b 2 ) 
log(QVA2) 



(7) 



We will take for the number of flavors nf = 5 and also Aqcp — 200 MeV. 
The correlation function $ a (x, k T ) is now defined as |p|,|l4| 



<f>(x,k T ) 



dk' 



(27T) 



(P, S\^(0) £+[0, -00] C + [-oo, x) ^(x)\P, S) 



(8) 



k+=XP+ Mr, 



with for each colored field a path-ordered exponential like 



C + [— 00, x] = V cxp —ig / dy A + (y) 



(9) 



The correlation function $> a (x, k T ) also has a factorization scale dependence (/x = bo/b), which is linked to that 
of the hard scattering part H. 

The correlation function $ a (x,fc T ) contains the nonperturbative dependence on the transverse momentum, 
which cannot be calculated perturbatively and should be fitted to experiment. In for example Refs. |^ |?J 
it is explained how these functions can be replaced by the ordinary parton distribution functions $ a (x), by 
introduction of a ^-regulator, e.g. the usual cut-off 6 max , via &» = b/y/l + 6 2 /6 2 nax , and by introducing a 
nonperturbative Sudakov factor Snp- Schematically, this proceeds as follows. The first term in Eq. (^) can be 
written as an overall b integration of an integrand called W(b), 



da 



d^ldx\dx2d 2 qT 



dx 



dx 



d 2 b 



,ib-q T 



W(b) + Y. 



(10) 



This function W(b) can be trivially rewritten as W(b*) x (W(b)/W(b*)). The first term (W(b*)) can be 
calculated within perturbation theory since a s (bo/b*) is always small (we will take the usual 6 max = 0.5 GeV -1 , 
hence a s (&o/fr* > 2) < 0.3) and the second term can be shown to be of the form 
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11 ih) = exp(-ln(Q 2 /Q 2 )g Q (b)-g A (x,b)-g B (x,b))=exp(-S NP (b)). (11) 



The functions g are not calculable in perturbation theory and need to be fitted to experiment. The functions 
gA(%, b) and gs(x, b) parameterize the intrinsic transverse momentum of the functions $> a (x, 6) and $ 6 (x, b) (in 
the polarized case <f>° is actually a function of b, see further comments below). It will be important later on 
that the functions g A and g p are independent of the scale Q. Finally, a lightcone expansion allows to express 
$ a (x, &*) in terms of $ a (x). 

We will consider $ a (x, b*)<I> b (x, b*)H{x, x; Q) only to lowest order in a s . The reason is that the perturbative 
tail of $°(x,fo*) is not known for the Collins function. Because of the unknown size of the distribution and/or 
fragmentation functions appearing in the asymmetries, the magnitude of the asymmetries cannot be estimated. 
Therefore, here we will be interested in the effect of soft gluons on the transverse momentum dependence of the 
asymmetries, which determines the Q 2 dependence of the Qt dependence. Although the hard scattering part 
is restricted to tree level, one can draw conclusions on the size of the effect of Sudakov factors by comparison 
to the tree level result, because both results will receive the same corrections to the hard part. Corrections to 
the hard scattering part can be easily inserted with hindsight (for a one-loop example cf. Ref. [[|, pages 432 
and 433). As said we are interested in the Q 2 evolution of the Qt dependence of the asymmetries. This Q 2 
dependence arises in a nontrivial way via Bessel functions, such that a numerical study is warranted. 

In Refs. [§-0 unpolarized scattering is considered, where only dependence on the length of b appears. Here we 
will consider polarized scattering and the above sketched procedure of Refs. |j-|7|] to factor a nonperturbative 
part from W(b), can also be applied after <E> (x,fc T ) is parameterized in terms of functions of k T . Explicit 
factors of k T -S T or fi T can be included in the hard scattering part H. 

The above factorization theorem for the Drell-Yan process provides the justification of the factorized ex- 
pression for the double transverse spin asymmetry studied in Ref. In Ref. ||^| an analogous factorization 
theorem had been discussed earlier for back-to-back jets in electron-positron annihilation, which comes down 
to a replacement of <&{x, k T ) — * A(z, k T ) -the fragmentation correlation function-. In Ref. it is argued 
that a similar factorization theorem holds for semi- inclusive DIS. Therefore, factorization theorems analogous 
to the one discussed above (differing only by obvious replacements), provide the theoretical foundation for the 
examples to be considered below. 

It is interesting to see what happens if one takes tree level everywhere in the above factorized expression Eq. 
(|J) (ignoring the Y term) . At tree level one finds the reduction 

/ iSf *- ih - {PT+hT - qT) e- S(b) - 5 2 (Pt + k T - q T ), (12) 
such that one arrives at the expression used by Ralston and Soper for the hadron tensor 15 16|JI|| 

W = \J d 2 p T d 2 k T S 2 (p T + k T -q T ) Tr($( 2 ; 1 ,p T )F/ I ¥(x 2 ,fc T )^)| p+ , fc - + (y^7)- (13) 

The vertex V/ 1 can be either the photon, Z or W boson vertex. In the above unpolarized case $(x,p T ) is 
parameterized as \ fi{x,p T ) f , where f\ is the parton momentum distribution function. Collins ]l7| has shown 
that polarization does not affect the factorization theorems and therefore, we conclude that the tree level 
formalism of Ralston and Soper Eq. ( |l3| ) is in accordance with the factorization theorem Eq. (Q) also in the case 



when polarization is taken into account. Clearly, Eq. (13) is only applicable for Qt values of the order of the 



intrinsic transverse momentum. To extend the formalism of Ralston and Soper to larger values of Qt, but still 
under the restriction of Q\ <C Q 2 , one can use the above factorized expression Eq. (|5|) (without the Y term) 
beyond tree level. 



III. UNPOLARIZED ASYMMETRY IN ELECTRON-POSITRON ANNIHILATION 

In this section we investigate a cos(20) asymmetry in the process of electron-positron annihilation into two 
almost back-to-back pions e + e~ — > ir + Tr~X. This process is very similar to the Drell-Yan process, so we 
can employ the same expression as given in Eq. (||), but with the correlation functions 4>, replaced by the 
fragmentation correlation function |18| 



4 



A«(*)=E / (^i e* z (0|^(x)/:-[x,oo]|P^;X)(P,5;X|/:-[^,0]^(0)|0), (14) 

which is used to define transverse momentum dependent fragmentation functions. 

The asymmetry requires observing the transverse momentum of the vector boson compared to the two pions. 
The tree level asymmetry expression for the cos(2</>) asymmetry in the process e + e~ — > 7r + 7r~ X was discussed 
in Refs. to which we refer for details. Here we shortly repeat the essentials. The asymmetry depends on 

the fragmentation function associated with the Collins effect [jl[ . The definition of the Collins function is 
given by ||) 

A(z,k T ) = JL L^k*) L +H i{z,k? T ) (15) 

where we have only displayed the fragmentation functions that can be present for unpolarizcd hadron produc- 
tion (Di is the ordinary unpolarized fragmentation function). The cos(20) asymmetry is an azimuthal spin 
asymmetry in the sense that the asymmetry arises due to the correlation of the transverse spin states of the 
quark-antiquark pair. On average the quark and antiquark will not be transversely polarized, but for each 
particular event the transverse spin need not be zero and moreover, the spin states of the quark and antiquark 
are exactly correlated. Subsequently, the directions of the produced pions are correlated due to the Collins 
effect and this correlation does not average out after summing over all quark polarization states. The Collins 
effect correlates the azimuthal angle of the transverse spin of a fragmenting quark with that of the transverse 
momentum of the produced hadron (both taken around the quark momentum), via a sin(</>) distribution of their 
difference angle <j). 

One finds for the leading order unpolarized cross section, taking into account both photon and Z-boson 
contributions EIT 



da(e + e- ->■ hih 2 X) 3a 2 



2^2 \ " 



dQd Zl dz 2 d 2 q T 02 Z 1 Z2 



IT 



Kf(y) T [£>i£>i] 



[K$(y) cos(2^) + Kl{y) sin(2^)] T 



(2h-p T h k T - p T -fc T ) ^yjjr 



(16) 
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= (.^ 2 + 




4 


= (.^ 2 - 
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The convolutions T will be dicussed below (Eqs. ( |35| ) and (]37|)), but first we will give the definition of the 
functions Kf (y) appearing in this expression and comment on the frame in which this cross section is expressed. 
As before, a is the flavor index; the functions Kf{y) arc defined as 

K?(y) = A(y) [e 2 a + 2g l v e a g vX i + c\c a lX2 \ - ^ [2g l A e a g a AXl + 4^X 2 ] , (17) 

mv) = B{y) [e 2 a + 2g v e a9vX i + c^Xa] , (18) 
Kt{y) = B(y) [2g v e a g AX3 ] , (19) 

which contain the combinations of the couplings 

(20) 
(21) 
(22) 

where gy and <?a are the vector and axial-vector couplings to the Z boson. The propagator factors are given by 

v - 1 Q 2 (Q 2 -m 2 z) m , 

Xl sin 2 (2^) (Q 2 ~ Ml) 2 + Y\M% ' 
l Q 2 

X2 = sm 2 (2e w ) Q 2 -M 2 z Xl ' (24) 
X3 = W^Ml^ (25) 
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The above cross section is expressed in the following frame. A normalized timelike vector i is defined by q 
(the vector boson momentum) and a normalized spacelike vector z is defined by P M = P M — (P ■ q/ 'q 2 ) q 11 for 
one of the outgoing momenta, for which we choose P2 (Pi are the hadron momenta), 



" Q ' 



P 2 -q 2 z 2 Q Q ' 



(26) 
(27) 



The azimuthal angles lie inside the plane orthogonal to t and z. In particular, <j> 1 gives the orientation of 



(g^ — iH u + z^z v ) l v , the perpendicular part of the lepton momentum I. 



In the cross sections we also encounter the following functions of y = l~/q~, which in the lepton center of 
mass frame equals y = (1 + cos#2)/2, where 9 2 is the angle of z with respect to the momentum of the incoming 
lepton I: 



i 



y + y' 



i 



1 + cos z 



0. 



B(y) = y (1 - y) c ^ -smH 2 , 



C(y) = (1 - 2y) 



COS 6*2 ■ 



(28) 
(29) 
(30) 



Since we have chosen P 2 to define the longitudinal direction, the momentum Pi can be used to express the 
directions orthogonal to i and z. One obtains P^ ± = P\ v (see Fig. |l|), 



Pf_L = -z\ q^ = z\ Q T h^, 
where we define the normalized vector = P?,/\P\±\. The angle 4>i is between h and l±. 



(31) 




FIG. 1. Kinematics of the annihilation process in the lepton center of mass frame for a back-to-back jet situation. 
P-2 is the momentum of a hadron in one jet, Pi is the momentum of a hadron in the opposite jet. 



The asymmetry A(q T ) at the Z mass (ignoring the photon and interference contributions, which can be easily 
included) is now defined by: 

da(e + e~ —>■ h\h 2 X) .„ , . . , ., 

d n d2ld z^ qT - {1 + ^^r)} , (32) 

with 

A confirmation of this asymmetry would confirm the Collins effect without the need of polarization. To avoid 
repeating irrelevant factors, we will first focus on the numerator and denominator of the term 



G 



K = 



Q% M 1 M 2 T [£>i£>i] 



(34) 



In Refs. tree level was considered, which means the expressions are valid in the region where the observed 
transverse momentum is very small compared to the hard scale (s), applicable only in the region of intrinsic 
transverse momentum. Eq. (p4) is a tree level expression if one uses the convolution notation 

T [DD] = J d 2 p T d 2 k T 8 2 {p T + k T ~ q T )D a (z 1 ,z 2 p 2 T )D a {z 2 , z\k 2 T ), (35) 

which only involves intrinsic transverse momenta. In order to go beyond this region, we include the Sudakov 
factor arising from resummed perturbative corrections to the transverse momentum distribution, by considering 
the approach of Ref. H|[. This will extend the range of applicability from the region of intrinsic transverse 
momentum to the region of moderate q T values (still under the restriction that q\ <C Q 2 )- This is outlined for 
the general case in the previous section; here we will apply it to the cos(2</>) asymmetry. 

Resummation of soft gluons into Sudakov form factors (see e.g. PQl) results in a replacement in Eq. (^5|) of 



S 2 (p T + k T -q T )^ J e -* (Pr+fcT - qr ) e -sW (36) 
leading to (suppressing the flavor indices) 

1 poo 

= — / dbbJ (bQ T )e- s ^D( Zl ,b)D{z 2 ,b). (37) 
27r Jo 

The function D denotes the Fourier transform of D. In order to compute the above expression, we will assume 
a Gaussian transverse momentum dependence for D±(z, z 2 k\\. 

D ± (z, z 2 k 2 T ) = D t (z) R 2 U exp(-^fe| ) /7r2 2 = D ^ qq^.rj/^ ( 38 ) 
For details see Ref. . Taking the Fourier transform of Eq. (^) yields 

D^b 2 ) = D t (x) exp M (39) 



The numerator in Eq. (p4[), 

d 2 b 



/d b 
(27r) 

x J d 2 p T d 2 k T (2q T -p T q T k T - q\p T -k T ) e ~^-(p T +k T ) D{z ll z\p 2 T )15{z 2 ,z 2 k 2 T ), (40) 



cannot be treated exactly like the denominator. A model for the transverse momentum dependence of the 
function is needed. In order to get a first estimate, or rather an upper bound for the asymmetry, one might 
be inclined to assume the maximally allowed function, by saturating the bound satisfied by pl[ 

\k T \ \H^{z, \k T \)\ < z M h D x {z, \k T \), (41) 

producing a l/|fc T | behavior of H±(z, k T ). However, this is not consistent with the fact that the Collins effect 
should vanish in the limit k T —> 0. 

A model by Collins jjj suggests the following transverse momentum dependence. Collins' parameterization 
for the fragmentation function is (note that Collins uses the function AD^/a ~ e^STikTj H\ > where St is 
the transverse spin of the fragmenting quark) 

D^z.k 2 ,) k\ + M 2 1 v ' v ;J z y ' 
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where Mh is the mass of the produced hadron and Mc is the quark mass that appears in a dressed fermion 
propagator i(A(k 2 ))t + B(k 2 )M c )/(k 2 - M£), the functions A and B are unity at k 2 = Mq. But for the 
present purpose, the additional fall-off with l/k T on top of the Gaussian fall-off is not needed. We will 
restrict to a Gaussian fall-off (with unknown magnitude) and assume the simple form for H^~ a (z, z 2 k' T ) = 
Hi a (z) Q{\k T \]R)j z 2 . Here one should take the radius R to be larger than R u of the unpolarized function D\, 
such as to satisfy the bound Eq. (|4l|) for all \ k T \. Also, we will assume that the fragmentation functions for both 
hadrons are Gaussians of equal width, i.e. we take R% = R2 = R and R u \ = R U 2 = R u and also Mi — M 2 — M. 
One finds 

J d 2 p T d 2 k T (2q T -p T q T k T - q 2 T p T -k T ) e" ,b ' (PT+ ^ ) Q(p T ;R) Q{k 2 T ]R) = 

i [2(9t • b)2 - 9 ^ 2] exp (- 2&) ' (43) 



which after application to Eq. (|40|) yields 

T[(2q T -p T q T -k T - q 2 T p T -k T ) DD] = J JA e ^(_J_) [ 2 (q T . b) 2 - q 2 T b 2 ] e' 3 ^ D( Zl ,b)t(z 2: b). 

(44) 

It is important to note that the factor 1/R 4 stems from the intrinsic transverse momentum of the functions 
Hi(z, k T ) and therefore, is not dependent on the scale Q (see discussion after Eq. (|i"i"|)). 

In analogy to A\ T , Eq. (HJ) can be transformed into (keeping in mind that in the end the summation over 
a, b in numerator and denominator should be performed separately) 

._ HfMHt-M A(QT)= ^y)i^y) A(QT% (45) 



4M'R' D*( 2l ) D»( 22 ) " D{( 2l )D!(«2) 

where for future reference we have also expressed the asymmetry in terms of the function (z) = 

z 2 J d 2 k T k T /(2M 2 ) H^{z, z 2 fe^) which occurs frequently in Qr-weighted cross sections (e.g. Q); furthermore, 

Mn v = M 2 r dbb 3 J 2 (bQ T ) exp (S(K) -S NP (b)) 

- ^dbbJ {bQ T ) exp(-S(b.)-S NP (b)) ' [ ' 

This is similar to the factor A(Qt) of Ref. Q], but with the replacement of Jo{bQr) — ► -hibQr) in the numerator. 
We note that unlike the case of A}f T investigated in Ref. Q, the present asymmetry does not need to oscillate 
as a function of Qt- Rather the vanishing of the numerator after integration over d 2 q T is due to the angular 
integration. Note that k only depends on Qt and not on fa. Also, we note that this asymmetry has a kinematic 
zero at Qt — 0, since h cannot be defined in that case (and indeed J 2 {Qt = 0) = 0). This is also seen from 
Ref. [p2[ , where a similar asymmetry factor has been investigated as a function of Qt- 

Here the main focus will be on the factor A(Qt), which is a measure for the effect of the Sudakov factors on 
the asymmetry compared to the tree level result exp [— (R 2 — R 2 )Q 2 ^/2] M 2 Q T R 6 /R 2 (cf. Eq. (|4g|)), which is 
valid only for values of Qt of the order of the intrinsic transverse momentum. 

In the above expression we have introduced the usual cut-off 6 max , via 6* = b/y/T+ fr 2 /^max an( i replaced 
\b 2 jR 2 — > Sjyp(b), for which we take one of the standard nonperturbative smearing functions, needed to 
describe the low qr region properly. It is important to realize that 5jvp(6) is introduced only in part to take 
care of the smearing due to the intrinsic transverse momentum (cf. Eq. (pi)), hence one cannot simply equate 
i b 2 /R 2 with Sivp(b). But taking into account the term exp (—5 b 2 /R 2 ) in addition to SNp(b) will just produce 
a change in the coefficient of the b 2 term in S]yp(b). To keep the unpolarized cross section unaffected, we will 
therefore introduce as nonperturbative term exp [— S^pib) + \ b 2 /R 2 ~\ , in order not to count the contribution 
from intrinsic transverse momentum twice. It is also worth mentioning again that since R 7^ R u (important at 
tree level), SNp(b) need not be the same in numerator and denominator (much less relevant however, since it 
affects only gA,gB and not gQ, cf. Eq. (|Tl|) ) and that in principle, it can depend on Z\ and z 2 , but we will not 
take into account these refinements. 



S 



Here we will take for the nonperturbative Sudakov factor the parameterization of Ladinsky-Yuan (Ref. ES 
with X1X2 — 10~ 2 ), 



SNp(b) = gib 2 + g 2 b 2 In 



_Q_ 

2Qo 



(47) 



with g x = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , Q = 1.6 GeV and 6 max = 0.5 GeV" 1 . Part of the results will depend 
considerably on this choice and this issue will be addressed in detail below. The value g\ = 0.11 GeV 2 can be 
viewed as an intrinsic transverse momentum of (p T ) = l/R 2 = (220 MeV) 2 . 

The reason we have chosen the parameterization of Ladinsky-Yuan |2^| , which is fitted to the transverse mo- 
mentum distribution of W/Z production in pp (pp) scattering, is that unfortunately there is no nonperturbative 
Sudakov factor available for e + e _ — > A + B + X, except for the energy-energy correlation function obtain from 
low energy data p4| ]. This is surprising considering the wealth of data from LEP experiments. Moreover, the 
nonperturbative Sudakov factor in the energy-energy correlation function (cx J2a b I dz\Z\dz-iZiQ 2 do j dQ\) 
as fitted in Ref. is not the same as the one in the differential cross section as a function of the lightcone 
momentum fractions. For related discussions, also relevant for SIDIS, see Refs. p5| , p6| . 

Although the nonperturbative Sudakov factors for Drell-Yan and e + e~ need not be related, the Snp of Ref. 
p3| can be viewed as a generic one and allows us to study the general features of the Sudakov suppression. 
We will also investigate the dependence on the nonperturbative Sudakov factor by varying the parameters. 
At a later stage one can always insert a more appropriate (phcnomcnologically determined) nonperturbative 
Sudakov factor into the asymmetry expression Eq. ([46]). Our results also underline the importance of a good 
determination of Snp- 

In Fig. U the asymmetry factor A(Qt) is given at the scales Q = 30 GeV, Q — 60 GeV and Q = Mz. 



A 




J 1 2 3 4 5 6 

Q T [GeV] 

FIG. 2. The asymmetry factor A(Qt) (in units of M 2 ) at Q = 30 GeV (upper curve), Q — 60 GeV (middle curve) and 
at Q = 90 GeV. 



The maximum of the asymmetry factor A{Q T ) at Q = 30 GeV (S NP (b) = 1.41 fo 2 ), Q = 60 GeV (S NP (b) = 
1Mb 2 ) and Q = M z (S NP (b) = 2.05b 2 ) is seen to be 0.57 (at Q T - 3.6 GeV), 0.31 (at Q T - 3.8 GeV) and 
0.22 (at Qt ~ 4 GeV), respectively. One observes that the magnitude of the asymmetry factor goes down with 
increasing energy and the position of the maximum -and also the average Qt~ shifts to higher values of Qt- 

Now we will discuss the dependence of these results on the choice of 6 max and Snp- Taking a higher value 
of &max increases the Gaussian width. The above choice of 6 max = 0.5 GeV^ 1 can be considered as optimistic 
(1/frmax is the scale down to which one trusts perturbation theory), enhancing the asymmetry factor somewhat. 

The asymmetry factor decreases with increasing Gaussian smearing width in Snp, to which it has a consid- 
erable sensitivity. Empirically, we find that if the Gaussian width is reduced by a factor a, then the maximum 
of the asymmetry increases roughly by a and the corresponding value of Qt by y/a. In Fig. [5] this is illustrated 
for the asymmetry factor A(Qt) at Q — Mz- The solid curve is for SnpO 3 ) = 2.05 b 2 and the dashed curve for 
Snp(P) — 1.37 b 2 , where the latter width is taken from a recent two parameter fit |2ffl. 
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FIG. 3. The asymmetry factor A(Qt) (in units of M 2 ) at Q 
(b) = 1.37 b 2 (dashed curve). 



Q T [GeV] 

90GeV with S N p{b) 



2.05 b 2 (solid curve) and 



The decrease with energy of the maximum of the asymmetry was found not to be very sensitive to changes 
to the Gaussian width in Snp- We find that the decrease goes as Q~ 9 — Q~ l a . 

So far we have only considered the asymmetry as a function of Qt, which means the cross section needs to be 
kept differential in the angle <pi and the magnitude Qt of the transverse momentum q T . Since the asymmetry 
does not vanish after integration over Q\ and since the Q^-integrated cross section has been studied using 
LEP data ply , we will now consider that case. Note that the integration has to be done in numerator and 
denominator separately. Since the denominator has no dependence on <f>i (the cos (201 ) dependence belongs to 
the numerator) one can in fact integrate over q T completely. This corresponds to the inclusive cross section 
in which soft gluon contributions cancel, hence this integration over the denominator of the asymmetry factor 
A(Qt) yields 1 and has no dependence on Sudakov factors. 

The Q^-integrated result should be compared to tree level, therefore we will first give the tree level expression. 
The tree level asymmetry is (neglecting the small Cg dependent term in the denominator) 



A(°\Q T ) = 



Q 2 T R 2 exp(-R 2 Q%./2) s i n 2 2 £ q c a 2 Ht a {z x ) Hi a (z 2 ) 



AM 2 R 2 U exp(-i? 2 Q 2 /2) 1 + cos 2 9 2 £ fc c\ D\{ Zl ) D\(z 2 ) 



(48) 



Note that if one integrates the cross section over Qff,, one should be careful to retain all dependences on Qt in 
both numerator and denominator separately. From the above expression we also infer that A(Qt) itself should 
be compared with the tree level quantity exp [— (R 2 — R 2 l )Q 2 ^/2] M 2 Q^R 6 /R 2 as mentioned before. 



Tree level: R„ 2 =l, R 2 = 3/2*R„ 2 



Sudakov result (xlO) 




Q T [GeV] 



FIG. 4. The asymmetry factor A(Qt) (in units of M 2 ) at Q — 90 GeV multiplied by a factor 10 (solid curve) and the 
tree level quantity (in units of M 2 ) using = 1 GeV -2 and R 2 /R^ = 3/2. 



In Fig. Q we have displayed the comparison of A(Qt) at Q = 90 GeV and the tree level quantity using the 
values R 2 L — 1 GeV -2 and R 2 jR\ = 3/2, which were chosen such as to minimize the magnitude. The value 
i? 2 = 1 GeV -2 can be regarded as too small already. We conclude that inclusion of Sudakov factors has the 
effect of suppressing the tree level result by at least an order of magnitude. This is important to keep in mind 
when making predictions of transverse momentum dependent azimuthal spin asymmetries based on tree level 
expressions. 
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For the particular case of a produced ir + and ir~ , we will make the following simplifying assumptions. We 
assume D"~* 7r (z) = Df^ v (z), D^ 1 * (z) = D"^*' (z) and neglect unfavored fragmentation functions like 
£jd-nr anc j s j m ji ar iy f or the Collins functions. As a consequence of these assumptions the fragmentation 
functions can be taken outside the flavor summation and appear as a square (at the average Zi, z 2 ). After Q\ 
integration one arrives at a cross section differential in the angle 4>\ of q T : 



do 



dfld(f>x 







ll + A 









n 2 



F(y) cos(2^i) 



(49) 



where 



sin 2 2 £ Q (^ 2 



;c(j/)) 



i 



9 a 2 ) 



COS 2 02 



E fc (g b v 2 



9 b A 2 ) 



(50) 



which is largest at 6 2 = 90°: F max w -0.5 (for a plot of the full factor F(0 2 ) see Ref. [§). 

For the prefactor A one finds at tree level = l/(2M 2 R 2 ) = {p 2 T )/(2M 2 ) and one should assume a 
typical intrinsic transverse momentum squared value, in the range of (p 2 ^) w (200 — 700 MeV) 2 . For pions this 
means A<® w 1 - 12. This can be compared to A^ = 6/tt of Ref. (28), which yields (p|) w (270 MeV) 2 (or 
i? 2 Ri 13GeV -2 ), which seems to be a reasonable value (remember that R 2 must be larger than i? 2 ). 

After including the Sudakov factors our numerical calculation yields A — 0.07 (S^pQ), Q = Mz) = 2.05 b 2 ), 
much smaller than the tree level values discussed above. The Gaussian width in Snp is not crucial for this 
conclusion; if one takes a much smaller width of for instance 1 GeV 2 , then one obtains A = 0.11. It has to be 
emphasized that the width used in SnpQ), Q) should increase with increasing energy. 

Our result shows that upon including Sudakov factors one retrieves parton model characteristics (also noted 
in Ref. ||), but with transverse momentum spreads that are significantly larger than would be expected from 
intrinsic transverse momentum (this is supported by the presently available parameterizations of Sn p in various 
processes, see e.g. [p5|). This means that a large tree level value for A^ (e.g. 6/71"), would lead to an extracted 
Collins function that is considerably smaller than if Sudakov factors would have been included. After the 
nonperturbative Sudakov factor has been obtained from LEP data, one can estimate exactly how much. For 
our choice of SNp(b,Q — Mz) = 2.05b 2 , the effect on the magnitude of the Collins function is an additional 
factor -y/6/(0.077r) » 5 compared to the extraction using the tree level expression. 

We conclude that the Sudakov factors produce a strong suppression compared to the tree level result. A tree 
level analysis applied at Q = Mz is expected to overestimate the asymmetry at least by an order of magnitude 
and therefore, it will underestimate the Collins function significantly. Hence, Sudakov factors should be taken 
into account when extracting the Collins function from LEP data or in general, from e + e~ data obtained at 
high values of y/s. 



IV. COLLINS ASYMMETRY IN SIDIS 



The Collins function originally was shown to lead to a single spin azimuthal asymmetry in semi-inclusive 
DIS 0] . This asymmetry has received much attention, since it would provide an additional way of accessing 
the transversity function hi. A preliminary determination of the Collins asymmetry from SMC data has been 
performed p9| , yielding an asymmetry for ir + production of 11% ± 6%. Also, the sin</> azimuthal asymmetry 
as recently measured by the HERMES Collaboration 1 30 1 might be related to the Collins asymmetry, providing 
further indication that the Collins function is nonzero. Measurements by HERMES and COMPASS are expected 
to provide more conclusive information on the Collins effect and its magnitude. Since these experiments are 
performed at different energies, it is important to know the Q 2 dependence of the asymmetry. Here we will 
investigate the Q 2 dependence of the transverse momentum (Qt) distribution of the asymmetry. 

From Ref. |3l| we extract the expression for this single spin asymmetry (extended to include contributions 
from Z-boson exchange; the expressions for VT-boson exchange can be found in Ref. pl|): 



da(W -> t'hX) _ 2a 2 xz 2 s 
dxdzdyd<f>id 2 q T Q 4 



£ \KUy)^[.hD 1 }-\S T \KUy)sm( 



(51) 
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where S T is the transverse spin of the incoming hadron H . The couplings K^(y) and K^(y) are of the same 
form as before, except that now 



(52) 
(53) 
(54) 



My) = \ i - v + 2* 

B{y) = -{l-y), 
C(y) = -y(2-y) , 



and y = (P ■ q)/(P ■ I) ~ q /I (I is the momentum of the beam lepton). Also, since Q 2 is now space- like, the 
width r z can be ignored: 



1 



Xx 



Q 2 



sin 2 (20 w ) Q 2 



X2 = (xiY 



(55) 
(56) 



The azimuthal angle 4> e h (4>g) around the three-momentum of the virtual boson is between the lepton's three- 
momentum and the outgoing hadron's three-momentum (transverse spin), cf. Ref. pl[ |. 
If we write the cross section as 



da(£H -> I'hX) 
dxdzdydtfiid 2 q T 

the asymmetry analyzing power is given by 

A(q T ) 



{1-\S 2 



sin 



b{+4)A{q T )} 



(57) 



Y. a Kt{y)T[q T -k T h l Hi] 



1T ' Q T M h Y, a Kt{y)T[hD M 
Assuming again Gaussian transverse momentum dependence and including Sudakov factors one arrives at 



(58) 



A(q T ) 



E a K$(y) ht(x)Ht a {z) 
2M 2 R 2 j: b K\(y) fi{x)D\(z) 



A(Qt), 



where R 2 is the Gaussian width of and the asymmetry factor A(Qt) is dehned as 



A(Qt) = M h 



/ °° dbb 2 JiffiQr) ex P (S(K) -S NP (b)) 
J™dbbJ Q (bQ T ) exp(-S(b*)-S NP (b)) ' 



(59) 



(60) 



This should be compared with Eq. (|46| ) of the previous section. In Fig. || the asymmetry factor A(Qt) is given 
at the scales Q = 30 GeV, Q = 60 GeV and Q=M Z . 




FIG. 5. The asymmetry factor A(Qt) (in units of Mh) at 
at Q = 90 GeV. 



5 6 7 
Q T [GeV] 

= 30 GeV (upper curve), Q — 60 GeV (middle curve) and 
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The maximum of the asymmetry factor A(Qt) at Q = 30 GeV, Q = 60 GeV and Q = Mz is seen to be 0.60 
(at Q T ~ 3.1 GeV), 0.42 (at Q T ~ 3.4 GeV) and 0.34 (at Q T ~ 3.5 GeV), respectively. Again we note that the 
magnitude of the asymmetry factor goes down with increasing energy and the position of the maximum shifts 
to higher values of Qt- 

As before we have used the nonperturbative Sudakov factor by Ladinsky and Yuan, Eq. (^), in the absence 



of a well-established Snp for SIDIS (cf. Ref. [|25 26 1). But the results obtained here can nevertheless be viewed 
as generic, only changes in the specific numbers are expected. Like for the previous cos(2</>) asymmetry, the 
asymmetry factor decreases with increasing Gaussian smearing width. Empirically we find that if the Gaussian 
width in Snp is reduced by a factor a, then both the maximum of the asymmetry factor and the corresponding 
value of Qt increase roughly by ^fa. However, the decrease with energy of the maximum of the asymmetry is 
not very sensitive to changes in this Gaussian width. We find that the decrease goes as Q~ ' 5 — Q~ ' 6 , which is 
considerably slower than for the cos (20) asymmetry. In the next section we will discuss the comparison between 
different asymmetries in more detail. 

We will now compare these results to the tree level expression for this Collins asymmetry. At tree level the 
convolution T (where w denotes a weight function) is given by 

F[w(p T ,k T )fD] = J d 2 p T d 2 k T 5 2 (p T + q T - k T )w (p T ,k T ) f a (x,p T )D a (z,z 2 k T ) . (61) 

If Gaussian transverse momentum distribution and fragmentation functions are assumed, one obtains at tree 
level 

A (o )(n x Q T R 2 eM-R 2 Q T m Ea my) K{x)Ht a {z) 

Wt) 2M h RleM-R 2 uQ T /2) H b Kt{ V ) f\{x)D\{z) ■ 1 > 

Therefore, one should compare A(Qt) with the tree level quantity exp [— (R 2 — R 2 l )Q' T /2\ MhQrR 4 / Rt- 




Qt [GeV] 

FIG. 6. The asymmetry factor A(Qt) (in units of Mh) at Q = 90 GeV multiplied by a factor 5 (solid curve) and the 
tree level quantity (in units of M h ) using R-t = 1 GeV -2 and R 2 /Rl = 3/2. 



In Fig. ^ we have displayed the comparison of A(Qt) at Q — 90 GeV and the tree level quantity using 
again the values R^ = 1 GeV -2 and R 2 / R 2 t — 3/2, which minimize the maximum asymmetry value. As before 
we conclude that inclusion of Sudakov factors has the effect of suppressing the tree level result. It is also 
clear that for the Collins asymmetry (which arises from a single Collins effect) the Sudakov suppression is less 
severe than for the case of the cos(2</>) asymmetry (which depends on the Collins effect squared). One also 
observes a pronounced increase of the average Qt upon inclusion of Sudakov factors. The difference to tree 
level becomes even more pronounced as the choice of R\ is increased to larger values more appropriate for a 
tree level analysis. Hence, also for transverse momentum dependent azimuthal spin asymmetries in SIDIS one 
needs to include Sudakov factors in order to extract distribution and fragmentation functions reliably. 

We now turn to a more general discussion of the different types of transverse momentum dependent az- 
imuthal spin asymmetries, categorized by different transverse momentum weights and angular dependences. 
The comparisons are most cleanly done for the asymptotic Q 2 case. 
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V. ASYMPTOTIC BEHAVIOR 



In the previous sections we have studied asymmetries with different types of weights, namely 2 h p T h-k T — 
p T -k T and hk T . The former (and also p T k T in the asymmetry A^ T of Ref. B) is typical of a double 
transverse spin asymmetry. It has two powers of transverse momentum and is therefore expected to decrease 
faster than the (single spin) asymmetries which have weights with a single power. The reason for this difference 
in decrease is that the asymmetries, which are convolutions of transverse momentum distributions, are largest if 
these distributions are large and have a large overlap. The distributions are largest at small p T and k T and the 
overlap is largest at small transverse momentum qx- However, the explicit powers of transverse momentum in 
the weights tend to suppress the region of small transverse momentum. More powers of transverse momentum 
in the weight thus implies more suppression. The power n of b n in the numerator compared to the denominator 
is a measure of this effect. A weight proportional to p T k 3 T will lead to n — 2 (cf. Eq. ((46|)), whereas p l T or k 3 T will 
lead to n = 1 (cf. Eq. d6C])). Larger n means, larger suppression and since the transverse momentum distributions 
broaden and decrease in magnitude with increasing energy, the convolution also decreases. Asymmetries with 
n = 1 will fall off slower with energy and the Sudakov suppression is less severe. This natural expectation is 
clearly observed in the numerical analysis of the examples investigated here and in Ref. [Q . 

It is interesting to see what is the asymptotic behavior of the Sudakov suppression. As is usually done, we 
will assume that at high values of Q the nonperturbative factor is irrelevant (true only asymptotically) and 
consider asymmetry factors of the form 

- 4 "< m(QT) = M r^J ( & Q T )exp(-^)) ■ (63) 

The asymptotic behavior of the denominator at Qt = (relevant for the unpolarizcd differential cross section) 
has been studied before J32|j7| by means of a saddle point approximation. Unfortunately, such a saddle point 
approximation cannot be applied to the ratios A nim (QT) for m ^ and Qt 0, since one cannot take 
Jm(bQr) ~ Jm(bspQT) to deduce the Q 2 dependence. For instance, since bspQr is in general not small, i.e. 
<C m (even under the assumption Q T <C Q 2 ), no approximation to the Bessel functions can be used in the 
region of interest (i.e., around the value of Qt where the asymmetry is maximal). 

In Ref. the asymmetry factor ^2,0 (Qt) arises, which does allow for a straightforward saddle point ap- 
proximation at Qt = 0, where the asymmetry is maximal. Therefore, we simply assume that the weight 
2h p T h k T — p T k T has similar behavior as p T k T . From our numerical studies we conclude that the maxi- 
mum of ^2,2(Qt) actually decreases faster with energy than the maximum of ^2,o(Qt) (at Qt = 0). Moreover, 
for the asymptotic behavior of -4 n ,o(QT = 0) one can obtain an analytic expression. First we would like to 
comment on the saddle point approximation of expressions of the form 

£V 6* exp (-£(&)) = J ^^lexp{(l + |)ln(6 2 A 2 )-5(fe)}. (64) 

As before we will keep only the leading term in the expansion of A(a s (/j,)) and a s (/i)/7r = l/(/?i ln(/i 2 /A 2 )): 

S(b,Q)= [ Q ^—^—-\n%. (65) 
J b 2 /b2 M 2 /3i lnO^/A 2 ) y 2 

The integral Eq. ( |64| ) has a saddle point at 

= X I A ) ■ (66) 



For n = we retrieve the power —0.41 of Refs. |32j,|7|. For general n one finds in the saddle point approximation 
(defining 7 „ = -C F /[C F + (1 + §)/?i]) 

d6 2 6«exp(-5(6)) S o< P ^J . (67) 

In this way one finds for the denominator of A2,o (Qt — 0) an approximate power behavior of (Q 2 )~ - 94 and for 
the numerator (Q 2 )~ 0,62 , when one takes Cf = 4/3 and /3i = 23/12 (corresponding to 5 flavors, but 6 flavors 
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only makes a few percent difference). The approximated asymmetry factor A2,o(Qt = OVthen has an asymptotic 
power behavior of (Q 2 ) -0 ' 32 , which is a similar decrease as found numerically in Ref. Q and apparently holds 
also for very large values of Q. At lower energies, the saddle point is not very pronounced and the introduction 
of the Gaussian smearing and 6* will have an effect on the approximation, hence we view this agreement as a 
coincidence. Of course, the saddle point approximation should become better with increasing Q 2 . 

For the ratios A n ,a(QT — 0) one can also derive an analytic expression for the asymptotic Q 2 dependence: 



f °° db 2 b n exp (S(b)) _ (M 2 b 2 y (Q 



Cp 

2\ TT 



In c n 



where c„ = (1 + C F /j3i)/(l + n/2 + C F /Pi)- For n = 2 this yields (Q 2 )" ' 32 for 5 flavors, which shows that 
the saddle point approximation is an extremely good approximation (the power of Q 2 differs by 0(1O~ 9 )). One 
now expects that the maximum value of A\.\{Qt) (appearing in the Collins asymmetry) to fall off at least as 
fast as Ai^{Qt — 0) (for all Q), which asymptotically goes as (Q 2 )~° ls . We note that the introduction of the 
Gaussian smearing and 6* have a considerable effect on the approximation at lower energies as our numerical 
studies demonstrate (in the studied range of 10 — 100 GeV we found approximately l/Q for max(A2,2(QT)) 
and 1/vQ for inax(Ai a (Qt)))- Our main conclusion is that for very high Q transverse momentum dependent 
azimuthal spin asymmetries will fall off as a fractional power of 1 /Q and that the behavior (or an upper bound 
to it) can approximately be found by looking at the power of b in the integrand. 

Finally, we want to mention that the increase of the Sudakov suppression with energy provides a solution to 
the following problem. In general, azimuthal spin asymmetries in collinear configurations, where partons are 
collinear to the parent hadrons, are suppressed by explicit powers of the hard energy scales. The transverse 
momentum factors in the weights have to be generated in the hard scattering part and for this one pays a price 
in terms of inverse powers of the hard scale. The collinear case can be viewed as averaged over the transverse 
momenta, relevant for qr T -integrated cross sections, which receives its main contribution from the region of 
small transverse momentum. On the other hand, in the case of q T -dependent differential cross sections at 
small q T (compared to Q), small parton transverse momentum needs to be included (one deals with nearly 
collinear partons therefore). In this case there need not be an explicit power suppression to generate similar 
azimuthal spin asymmetries (the average transverse momentum is now a scale in the problem, so dimcnsionless 
ratios using Q F rather than Q can be formed). This seems counterintuitive, since by considering a differential 
cross section for a nearly collinear configuration, one can obtain information without power suppression, that 
is power suppressed in the g T -integrated cross section, which receives its main contribution from the region of 



small transverse momentum also (for an example, see Ref. 33|). But we find that instead of explicit power 
suppression the presence of Sudakov factors gives rise to partial power suppression. In this sense the power 
suppression is effectively replaced by a Sudakov suppression. This means that the q T -dependent azimuthal 
spin asymmetries vanish with increasing energy, as do their q T -integrated counterparts. Hence, departing from 
the collinear configuration (that is, by including small transverse momenta) does not allow one to measure 
azimuthal spin asymmetries in (differential) cross sections at arbitrary high energies. 



VI. CONCLUSIONS 



In this article we have shown by quantitative examples how transverse momentum dependent azimuthal 
spin asymmetries are suppressed by Sudakov factors, in the region where the transverse momentum is much 
smaller than the large energy scale Q 2 . Physically, the Sudakov suppression stems from broadening of the 
transverse momentum distribution due to recoil from soft gluon radiation and the suppression increases with 
energy. This implies that tree level estimates of transverse momentum dependent azimuthal spin asymmetries 
tend to overestimate the magnitudes and increasingly so with rising energy. 

The size and Q 2 dependence of the Sudakov suppression have been studied numerically for two such asym- 
metries, both arising due to the Collins effect. The size of the suppression (compared to tree level) depends 
considerably on the nonperturbative Sudakov factor that must be determined from experiment; however, the 
Q 2 dependence of the suppression turned out to be much less sensitive to the nonperturbative input. 

We observe that in general the larger the power n of the transverse momentum in the weight of an asymmetry, 
the larger the suppression. For the Collins effect driven cos(20) asymmetry in electron-positron annihilation into 
two almost back-to-back pions (n = 2), the Sudakov suppression was numerically found to be approximately 
l/Q for the maximum of the asymmetry and an upper bound for the asymptotic behavior was found to be 
l/Q ' 6 . For the Collins effect single spin asymmetry in semi-inclusive deep inelastic scattering (n — 1), the 
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Sudakov suppression was numerically found to be approximately 1 / \[Q for the maximum of the asymmetry and 
an upper bound for the asymptotic behavior was found to be 1/Q 0A . For the maximum of the asymmetries 
of the type ^4 n ,o(Qr), the asymptotic Q 2 dependence could be calculated analytically. This provides upper 
bounds on the fall-off of A n , m {QT) > with m^O. 

Since our results depend on the input for the nonpcrturbative Sudakov factors Snp, which is not (well) deter- 
mined for the above processes, the numerical conclusions about the size and Q 2 dependence of the suppression 
should be viewed as generic, not as precise predictions. Therefore, we would like to stress the need for an 
extraction of the nonperturbative Sudakov factor from the process e + e~ — > A + B + X (for any two, almost 
back-to-back hadrons A and B) and from SIDIS. Considering the wealth of data from the LEP and HERA 
experiments this should pose no problem. 

The Sudakov suppression of the transverse momentum distribution of azimuthal spin asymmetries is significant 
already for Q values in the range of 10 — 100 GeV as can be concluded from comparison to tree level. In single 
spin asymmetries the magnitude of the suppression is less severe than in double spin asymmetries and in both 
cases a pronounced shift of the average Qt to higher values is observed. We conclude that it is essential to take 
into account Sudakov factors in transverse momentum dependent azimuthal spin asymmetries. 
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